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ABSTRACT 

Black-hole masses of the highest redshift quasars (4 < z < 6) are estimated using a previously presented 
scaling relationship, derived from reverberation mapping of nearby quasars, and compared to quasars at 
lower redshift. It is shown that the central black holes in luminous z > 4 quasars are very massive (> 
10 9 Mq). It is argued that the mass estimates of the high-z quasars are not subject to larger uncertainties 
than those for nearby quasars. Specifically, the large masses are not overestimates and the lack of similarly 
large black-hole masses in the nearby Universe does not rule out their existence at high-z. However, AGN 
host galaxies do not typically appear fully formed and/or evolved at these early epochs. This supports 
scenarios in which black holes build up mass very fast in a radiatively inefficient (or obscured) phase 
relative to the stars in their galaxies. Additionally, upper envelopes of A/bh ~ 10 10 Mq and Lboi ~ 10 48 
ergs s" 
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1. INTRODUCTION 

High-redshift quasars are of immediate interest to both 
cosmologists and active galactic nuclei (AGNs) specialists 
for several reasons. For example, quasars can be observed 
out to greater distances and with more ease than quiescent 
galaxies, owing to their unparalleled energy output over es- 
sentially the entire electromagnetic spectrum (e.g., Elvis 
et al. 1994). Since quasars reside in the central regions of 
massive galaxies (e.g., Boroson & Oke 1982; Boroson, Oke, 
& Green 1982; Kukula et al. 2001), they provide a valu- 
able beacon to locate galaxies at early times, which is diffi- 
cult by other means 1 (Lanzetta et al. 1999). Furthermore, 
quasar activity is intimately connected with the hierarchi- 
cal buildup of galaxies (e.g., Haehnelt, Natarajan, & Rees 
1998; Cattaneo, Haehnelt, & Rees 1999; Ferrarese & Mcr- 
ritt 2000; Gebhardt et al. 2000a; Haehnelt & Kauffmann 
2000; Colina et al. 2001; Menou, Haiman, & Narayanan 
2001). In fact, quasars possibly play a profound role for 
the formation and evolution of most, if not all, galaxies 
(e.g., Kauffmann & Haehnelt 2000; Milosavljevic & Mer- 
ritt 2001, 2002; Benson et al. 2002). The quasar phase 
may even precede the formation of the galaxies themselves 
(Rix et al. 2001; Kauffmann & Haehnelt 2000; Omont 
et al. 2001) and impact their early formation and evolu- 
tion (e.g., Silk & Rees 1998; Haiman, Madau, & Loeb 
1999; Theuns, Mo, & Schaye 2001; Di Matteo et al. 2003). 

By virtue of being associated with the highest density 
peaks in the matter distribution (e.g., Efstathiou & Rees 
1988; Volonteri, Haardt, & Madau 2002), quasars are an 
easy probe of these peaks at the earliest epochs and hence 
of some of the first structures formed in the early Universe 
(e.g., Haiman & Loeb 1998). In addition, quasars help 
probe the nature of the medium in which they are em- 



bedded. For example, the spectra of z rj 6 quasars place 
constraints on the epoch of re-ionization (Fan et al. 2002). 
The presence of luminous, massive 2 quasars at redshift 6, 
when the Universe was less than 10% of its present age, 
already place much-needed constraints on models of the 
earliest structure formation (Turner 1991; Haiman & Loeb 
2001; Volonteri et al. 2002). 

High-redshift quasars are also studied (a) to track the 
growth and evolution of the central black holes (Haiman 
& Loeb 2001; Merritt & Ferrarese 2001a; Ferrarese 2002; 
A. Steed et al. 2003, in preparation), (b) to under- 
stand the central engine (e.g., Peterson, Polidan, & Pogge 
2001; Peterson 2002) responsible for their monstrous en- 
ergy release, and (c) to understand the evolution of the 
quasars and active galaxies in general (e.g., Pei 1995; Boyle 
et al. 2000; see Hartwick & Schade 1990 and Osmer 2003 
for reviews). 

It is well established that quasars and their activity 
reign 3 between redshift 2 and 3 (Osmer 1982; Hartwick 
& Schade 1990; Warren, Hewett, & Osmer 1994; Schmidt, 
Schneider, & Gunn 1995; Kennefick, Djorgovski, & Mey- 
lan 1996; Fan et al. 2001a). The decline in comoving space 
density above redshifts of ~ 3.5 is significant not only 
for the optically selected quasars but also for the radio- 
selected subset (Shaver et al. 1996; Hook et al. 2002) and 
for X-ray selected sources (Hasinger 2002); the latter may 
decline more slowly with increasing redshift, but the statis- 
tics are still relatively poor. The question naturally arises 
whether or not the space density drop toward higher z 
also signifies a change in cardinal properties of the quasar 
population at that epoch, given the relatively short time 
available to generate massive black holes; at z > 3 the 
Universe was only <10% of its present age for a A = 
Universe (and <15% of its present age for A = 0.7). 



1 However, SI RTF, JWST, FIRST/ Herschel and ALMA will undoubtedly change this (e.g., Neufeld 1999; Omont 2003). 

2 This is typically estimated via the Eddington argument but these high black-hole masses are also confirmed in this work. 

3 Quasars were most numerous at the epoch at z m 2.5, yet they are still rare objects (e.g., Hartwick & Schade 1990) relative to the galaxy 
population (e.g., Huchra & Burg 1992). 
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Recent developments allow virial black-hole masses in 
distant quasars to be estimated using scaling relationships 
calibrated to reverberation mapping masses of nearby 
AGNs and quasars (Wandel, Peterson, & Malkan 1999; 
Kaspi et al. 2000; Vestergaard 2002, hereafter Paper I). 
In this work, these scaling relationships are used to es- 
timate the black-hole masses in high-z quasars. The 
main focus in this work is to compare the typical cen- 
tral masses of quasars just below and above the epoch 
z w 3.5 at which the comoving space density of quasars 
declines. In addition, these masses are discussed in light 
of our current knowledge of the nature of their host galax- 
ies. Representative quasar samples in the redshift range 
1.5 < z < 6.3 are studied (§ 3), including part of the 
large, high-quality data bases of high-redshift quasars now 
available. Several surveys have in the past ~10 years 
turned up a large number of previously unknown z w 4 
quasars (e.g., Kcnncfick ct al. 1995a, 1995b, 1996; Storrie- 
Lombardi et al. 1996; Schneider, Schmidt, & Gunn 1997; 
see also Constantin et al. 2002). The Sloan Digital Sky 
Survey (SDSS) has in recent years had particular success 
in uncovering an even larger number of quasars at redshifts 
greater than 3.5 (Fan et al. 1999, 2000, 2001b, 2001c; An- 
derson et al. 2001). For a large subset of these quasars for 
which spectra are published, the black-hole masses were 
estimated in this work and compared with mass estimates 
of lower-redshift quasars, namely the well-selected z « 2 
quasar sample of Vestergaard (2000). Given current and 
earlier survey capabilities, these samples necessarily con- 
sist of luminous quasars. Since data on nearby luminous 
quasars and AGNs readily exist in the Bright Quasar Sur- 
vey sample (hereafter BQS; Schmidt & Green 1983), this 
sample was also included for comparison and reference. 
The lack of objects between 0.5 < z < 1.5 does not affect 
the conclusions of this study. 

The structure of this paper is as follows. The approach 
in estimating Mbh, -^boi, and -Lboi/^Edd is outlined next 
(§ 2). The data are described in § 3, presented in detail in 
§ 4, and discussed in § 5. In particular, since the central 
masses are very large (~ 10 9 M Q ) for the z > 4 quasars, 
the efficacy of the scaling relations is examined in § 5.1; it 
is concluded that the scaling relations are fully applicable 
to the high-z quasars with an uncertainty of a factor <4. 
Section 5.2 discusses the fact that although the masses of 
these high-z active black holes are rather large, the AGN 
host galaxies at z > 4 appear to be rather young, and 
even small, young galaxies at z « 3 are capable of host- 
ing massive black holes. This leads to the conclusion that 
super-massive black holes were in place long before the 
bulk of the stellar mass in the quasar host galaxies. A 
summary and conclusions are provided in § 6. 

A cosmology with H = 75 km s _1 Mpc _1 , q = 0.5, and 
A = is used throughout, although the main conclusions 
of this study do not depend on the specific choice of cos- 
mology. However, it is worth mentioning that adoption of 
a cosmology of H = 70 km s _1 Mpc _1 , Qm — 0.3, £1a = 
0.7 will increase the mass estimates at the highest redshifts 
by no more than a factor 2.4 and by a decreasing factor at 
lower redshifts. 



2. M B H, Lbol, AND Lbol/LEdd MEASUREMENTS 

The samples of quasars studied here cover the redshift 
ranges z < 0.5 and 1.5 < z < 6.3. Therefore, the observed 
optical spectra cover restframe wavelengths including H (3 
for the BQS and the region around Civ for the higher 
redshift quasars. 

The virial masses, M cx v 2 R, are estimated from single- 
epoch spectral measurements using the relationships pre- 
sented in Paper I (eqns. 1, 2, A5, and 7), since they are 
relevant for spectroscopic data covering the wavelength re- 
gions around the H (3 and C iv lines. Here v is the Dopplcr 
width of an emission line and R is the radius of the line- 
emitting region. Note that also the high-ionization line- 
emitting gas, including C iv, yields virial masses consistent 
with those derived from H/3 in all the objects for which 
this can be tested (Peterson & Wandel 2000; Onken & 
Peterson 2002). Hence, the Civ line width is valid as a 
measure of v (see discussion in § 5.1). 

The important scaling relationship for the mass equa- 
tions listed below is the radius — luminosity relation, 
that allows the size of the broad line region (BLR) to be 
estimated from the optical continuum luminosity of the 
source. The empirical scaling law, applicable to Balmer 
lines and based on reverberation-mapped nearby AGNs 
and quasars, is (Kaspi et al. 2000): 

flBLR = <32.9 ± , 5) r Ai > (510oAr ' ' : 



10 44 



ergs s~ 



It — days. 
(1) 

The mass equation based on optical data (hereafter, the 
'optical relationship') is 



log Mbh/M q = 6.7 + 



(2) 



log 



/FWHM(H/3)\ 2 /AL A (5100A) 
\1000 km s- 1 ) V10 44 ergs s" 1 



0.7" 



and the mass equation relevant for UV data (the 'UV re- 
lationship') is 



log M BH /M 
log 



6.2 + 



/FWHM(C IV) \ / Ai A (1350A) \ 
V 1000 km s- 1 J \10 44 ergs s" 1 ) 



Q.T 



(3) 
, (4) 



where the sample standard deviation of the weighted mean 
is 0.45 and shows the intrinsic scatter in the sample (Pa- 
per I). This is confirmed by the black- hole masses of rever- 
beration mapped AGNs which exhibit an intrinsic scatter 
of about 0.3 — 0.5 dex (Peterson 2002) around the Mbh — o" 
relationship, established by local quiescent galaxies (e.g., 
Ferrarese & Merritt 2000; Gebhardt et al. 2000a; see also 
Gebhardt et al. 2000b and Ferrarese et al. 2001 for com- 
parisons of some AGN masses with the Mbh — cr relation). 
Both the optical and the UV relationships were calibrated 
(Paper I) to virial mass estimates of low-redshift active 
galaxies determined from reverberation-mapping analysis. 

The velocity dispersion, v, of the broad line gas is es- 
timated using the broad emission line width as charac- 
terized by FWHM. While equation (1) is valid for the 
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Balmer lines only, a similar scaling law, -Rblr °c L 0,7 , 
is assumed for the UV relationship; the constant of pro- 
portionality is determined by the calibration and included 
in equation (4). Note that when applying these scaling 
relationships to higher- z quasars, a large extrapolation in 
luminosity is not performed, which might otherwise raise 
concern. The efficacy of the mass scaling relationships is 
discussed in § 5.1. 

The bolomctric luminosities are estimated based on ex- 
isting average spectral energy distributions (SEDs). SEDs 
are only available for a small group of distant quasars. 
Kuhn et al. (2001) found no significant differences between 
the restframc optical and UV energy distributions for 15 
z > 3 quasars and low-z quasar SEDs. It is fair to assume 
the Kuhn et al. sample is representative, also of quasars at 
higher rcdshifts yet, given the similarity of quasar spectra 
at all redshifts (see § 5). Specifically, this means that the 
empirical average SED of low- z active galaxies and quasars 
(Elvis et al. 1994) is also representative of the average SED 
of more distant quasars. Here, the average radio-quiet 
quasar (RQQ) SED of Elvis et al. (1994) is updated to re- 
flect the more representative optical-to- X-ray slope of 1.43 
determined by Elvis, Risaliti, & Zamorani (2002) using 
new and improved X-ray data. This updated RQQ SED 
was used to estimate bolometric correction factors applica- 
ble also to radio-loud quasars (RLQs), because (1) high- z 
RLQs appear to constitute a similar low fraction (~10%) 
of all quasars as seen for z < 2 quasars (Stern et al. 2000), 
and (2) the average RLQ SED has not yet been updated 
(M. Elvis, 2002, private communication). The bolometric 
correction factors to the monochromatic continuum lumi- 
nosities, \L\, at rest wavelengths of 1350A, 1450A, and 
4400A were computed 4 to be 4.62, 4.65, and 9.74, respec- 
tively. Applying a bolometric correction to the UV contin- 
uum luminosities for the high-z quasars reduces the uncer- 
tainty in the Lboi estimate compared to using a bolomet- 
ric correction to luminosities at other wavelengths. The 
reasons are that one eliminates the uncertainty associated 
with the measured UV continuum slope and extrapolation 
to other wavelengths, and since A1350 and A1450 are close 
to the peak of the big blue bump, the associated bolo- 
metric correction factors are the smallest relative to other 
wavelengths. 

The conservative measurement uncertainties estimated 
from the line widths, continuum fluxes, and continuum 
slopes were propagated to estimate the measurement er- 
rors on Mbh, -^boi) and the Eddington ratios I/boi/^Edd oc 
-^boi/^BH- These uncertainties do not include the uncer- 
tainty associated with the use of bolometric corrections or 
the use of the simple scaling relationships in eqn. (3), (1), 
and (4) to estimate the Mbh- A statistical lcr uncertainty 
of about 0.4 dex and 0.5 dex (i.e., a factor 2.5 and 3 in the 
mass) for the optical and UV relationships, respectively, 
was found for the sample studied in Paper I. This addi- 
tional uncertainty must be kept in mind, especially when 
considering the masses in absolute terms. However, since 



Mbh can likewise be reasonably estimated for high red- 
shift quasars (§ 5.1), statistical analysis and comparison 
of the relative masses of high-redshift quasars should be 
reasonably sound. 

3. DATA 

Several samples of quasars at low, intermediate, and 
high redshift arc analyzed here as described in detail be- 
low. 



3.1. Low Redshift Quasars 

This sample consists of the 87 BQS quasars at z < 
0.5. For the objects with available reverberation-mapping 
masses, these masses are used, as computed in Paper I. 
For the remaining sources, the masses are estimated from 
equation (3). The measurements of line widths and contin- 
uum luminosities used are explained in detail in Paper I. 
Briefly, the FWHM(H 0) measurements were primarily 
adopted from Boroson & Green (1992). The monochro- 
matic luminosities, L\ (5100A), were determined for most 
of the objects from the spectrophotometry by Neugebauer 
et al. (1987), but for a few sources £ A (5100A) was com- 
puted from 4400 A flux densities (Kellerman et al. 1989) 
by adopting a power-law slope, a = —0.5 (F„ oc v a ). 
Bolometric luminosities for different subsets of the BQS 
quasars have been measured by a few authors (Sanders 
et al. 1989; Elvis et al. 1994; Lonsdale, Smith, & Lons- 
dale 1995; Wilkes et al. 1999a). Comparison of the mea- 
surements for quasars common between two or more of 
these studies shows that the Lboi values are generally rel- 
atively similar 5 within a common scatter of 0.2 — 0.3 dex, 
consistent with the estimated measurement uncertainties 
quoted by Elvis et al. (1994). In the interest of simplicity 
and consistency with the data used for the quasar samples 
at higher-redshifts, the Lboi values were here estimated 
as (9.74 ± 4.3) x AL A (4400A) (§ 2). Adopting this scal- 
ing introduces an uncertainty similar to the measurement 
uncertainties quoted above and is thus considered accept- 
able. Errors on the parameters are propagated using the 
measurement uncertainties on the primary measurements 
as described in Paper I. 

3.2. Intermediate Redshift Quasars 

The sample of 114 quasars at intermediate redshifts 
(1.5 < z < 3.5) studied here is that analyzed by 
Vestergaard (2000). The sample consists of 68 RLQs 
and 48 RQQs. The RLQs were selected from Barthcl 
et al. (1988), Barthel, Tytler, & Thomson (1990), Murphy, 
Browne, & Perley (1993), and the 3C and 4C radio cat- 
alogs to span the range of observed quasar radio spectral 
indices; this allows a study of how the spectral features 
may change with source inclination with respect to our 
line of sight (e.g., Vestergaard, Wilkes, & Barthel 2000). 
Most of the RLQ data are from Barthel et al. (1990). The 
individual RQQs were selected from Hewitt & Burbidge 

4 The bolometric correction to AL^(4400A) of 9.74 is slightly lower than the average correction of 11.8 established by Elvis et al. (1994). The 
reason is in part the weaker contribution from the X-rays, and in part due to the use of the average SED in this work compared to the average of 
the individual bolometric correction factors quoted by Elvis et al. (1994). The scaling factor of 9.74 is preferred for AL^(4400A) for consistency 
with the 1350A and 1450A bolometric correction factors which were not quoted by Elvis et al. 

5 While one can argue that one of the studies has better X-ray data, then another has better IR data, and for individual objects, the SED of 
one study may be preferred over that of another. However, it is beyond the scope of this work to compile updated SEDs based on the best 
data now available. 
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(1993) to be well-matched and paired in redshift and lu- 
minosity to individual quasars in the radio-loud subset. 
The purpose was to study the rest UV spectral charac- 
teristics with radio properties (Vestergaard 2000). Given 
the careful sample selection, this sample of 116 quasars 
is here considered reasonably representative of the z w 2 
quasar population to be compared with the higher redshift 
quasars. The full sample, as opposed to the pair-matched 
subset, is studied here to improve the statistical signifi- 
cance of the comparison. 

The Civ line width, FWHM(Civ), was measured re- 
liably for this sample on a smooth representation of the 
emission profile, obtained by modeling with 2 to 3 Gaus- 
sian functions (following Laor et al. 1994), above a global 
power-law continuum fit to line-free windows (e.g., Fran- 
cis et al. 1991; Vestergaard & Wilkes 2001). The FWHM 
uncertainties were estimated by varying the continuum 
fit within the continuum flux rms and remeasuring the 
FWHM. The continuum fluxes were measured from the 
best fit power-law and the flux uncertainties were deter- 
mined as the rms scatter around this fit in the continuum 
fitting windows. The bolometric luminosity is estimated 
from L A (1350A), i.e., L bol = 4.62 xAL A (1350A) [§ 2]. 

3.3. High Redshift Quasars 

The high-z quasar sample consists of —150 z > 3.5 
quasars, mostly from the SDSS, selected from recently 
published work, as detailed below. For consistency, the 
spectra of these quasars were measured in a similar man- 
ner, as follows: FWHM (Civ) was measured relative to a 
reasonable, best-fit continuum determined using virtually 
line-free continuum windows. The continuum flux at rest 
frame wavelengths 1350A and 1450A were measured from 
this best-fit power-law continuum. Conservative measure- 
ment uncertainties on line widths, continuum fluxes, and 
spectral slopes were estimated based on the continuum rms 
noise, as described above for the intermediate- z quasars. 
The spectral slope measured here is only used when more 
accurate, published values are unavailable. Lboi was esti- 
mated as 4.65 xAL A (1450A) (see § 2). 

The source studies from which the individual quasars 
are drawn are listed below. In each case, the details of the 
analysis specific to the individual data set are outlined. 

Constantin et al. (2002): High-quality spectra of 
these z « 4 quasars were kindly provided by the authors. 
Constantin et al. corrected their data for Galactic red- 
dening using Schlegcl, Finkbeiner, & Davis (1998) Ay val- 
ues and the Cardelli, Clayton, & Mathis (1989) extinc- 
tion curve. Constantin et al. quote an uncertainty of 0.2 
dex in their AB(1450A)-based, monochromatic luminosi- 
ties. The L A (1450A) values measured here were recali- 
brated to match the AB-bascd L A (1450A) for the few ob- 
jects showing a difference larger than 0.2 dex between the 
two sets of luminosity measurements. The spectra of 40 of 
the 44 quasars could be measured. A few quasars are too 
strongly absorbed in the C iv emission profile for a reason- 
able measurement of the line width. One quasar exhibited 
no conspicuous emission lines other than Lya and Nv. 

Anderson et al. (2001): The digital SDSS spectra of 
these —125 quasars are publicly available. The best qual- 



ity spectra were selected for measurement. A few other 
poorer quality spectra were also measured. The spec- 
tra were first corrected for Galactic reddening using the 
Schlegel et al. (1998) E(B-V) values obtained from NED 
and the Cardelli et al. (1989) reddening curve. For a sub- 
set (-60) of the quasars Fan et al. (1999, 2000, 2001b) 
present AB(1450A) magnitudes to which the spectral mea- 
surements were recalibrated. 

Fan et al. (1999, 2000, 2001b): The best-quality 
spectra (kindly provided by X. Fan) of these SDSS quasars 
that include a full, measurable C IV line profile were se- 
lected for measurement. The continuum luminosities mea- 
sured here at 1350A and 1450A were recalibrated to the 
(dereddened) AB(1450A) magnitudes. Fan et al. define 3 
quasar subsamplcs. The set of 59 quasars analyzed here 
(hereafter the "Fan et al. sample") is drawn from these 
subsamplcs, as follows: (a) 34 of the 39 quasars that com- 
prise the well-defined "color-selected sample" (3.6 < z < 5; 
i* < 20 mag), complete within the survey area (see Fan 
et al. 2001b for details), (b) 12 quasars from the "bright 
subsample" (i* < 20.2 mag) that are not included in (a), 
and (c) 13 of the 18 "faint subsample" quasars (i* < 21 
mag). As mentioned, the selection is mostly based on 
data quality (but see below). Ten quasars were com- 
mon between the originally selected "Anderson et al. " 
and the "Fan et al. " samples. Since neither of these sam- 
ples are complete, each of those quasars only occur once 
in the "high-z sample"; each object was included in the 
sample which has the better spectrum. The data of the 
"color-selected quasars" are presented by Fan et al. (1999, 
2001b) and Schneider et al. (2001). The photometry of 
Fan et al. (2001b) and line widths measured in the high 
S/N spectra of Constantin et al. were used for two ob- 
jects common to these studies. Five quasars were excluded 
from the "color-selected sample": the spectra of three of 
these quasars did not include the C iv emission line. For 
another object (J020731. 68+010348.9) the poor spectrum 
quality precludes a FWHM(Civ) measurement. Finally, 
the spectrum of a fifth object, J225529.09-003433.4, could 
not be located in any of the references listed or in the 
SDSS archive. It will be shown in the next section that 
the "color-selected sample" is statistically similar in the 
distributions of Lboi, -^bh, and L^i/L^dd to those of the 
other high-z quasar samples (3.5 < z < 5.5), as expected 
if this is a representative subset of z s=a 4 quasars. 

Pentericci et al. (2002): The authors present near- 
IR spectra of two of the z « 6 quasars discovered in SDSS 
(Fan et al. 2001c), namely SDSSp J103027.10+052455.0 
(z = 6.28) and SDSSp J130608.26+035626.3 (z = 6.00). 
The C IV line widths are estimated from enlarged versions 
of the published spectral plots, while the luminosities are 
adopted from Fan et al. (2001c). The width of Civ mea- 
sured here for SDSSp J103027. 10+052455.0 is consistent 
with that quoted by Pentericci et al. within the errors. 
The C iv profile is truncated in the blue-most part of the 
wing. The FWHM estimate (=7000t™° kms" 1 ) used 
here for Mbh was obtained by measuring a reasonable ex- 
tension of the observed blue profile to the continuum level 
indicated. The (large) FWHM measurement uncertainty 
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log (M BH / MJ log (L B0L / erg S"») log (L B0L / L Edd ) 



Fig. 1. — Distributions of estimated (a) black-hole mass, Mbh, {b) bolometric luminosity, Lboli and (c) Eddington luminosity ratio, 
^bol/^Edd f° r different redshift bins. The top panel consists of Bright Quasar Survey objects only. The middle panel contains the intermediate- 
z sample; the radio-quiet quasars (RQQs) are shown shaded. The bottom panel unshaded histogram contains all the (z > 3.5) quasars in the 
samples described in § 3.3. The distribution of the faint subsample of Fan et al. (2001b) is high-lighted with diagonal shading. The histogram 
for the "color-selected sample" of Fan et al. (2001b) is shown with a light horizontal shading for comparison. 



adopted is based on a conservative estimate of the contin- 
uum uncertainty. 

4. DISTRIBUTIONS OF MASS, LUMINOSITY, AND 
EDDINGTON RATIO 

4.1. Number Distributions 

The histograms of masses, bolometric luminosities, and 
Eddington ratios are shown in Figure 1 for different red- 
shift bins. Figure 2 shows the distributions of emission-line 
widths with redshift and continuum luminosities. 

It is evident from Figure 1 that the intermediate and 
high-rcdshift (z > 1.5) quasars are all massive (> 10 8 M Q ) 
and luminous (> 10 46 ergss~ 1 ) and have somewhat simi- 
lar Mbh and Lboi distributions. As discussed in § 4.2 the 



lower boundaries are due to selection effects. The two ra- 
dio types of the intermediate-redshift quasars, which are 
matched in luminosity and redshift, display no significant 
Mbh difference, in contrast to some studies claiming black- 
holes of RLQs are more massive (e.g., Laor 2000; McLure 
& Dunlop 2001). The radio-loud z ss 2 sample does have 
a larger number of luminous quasars than the radio-quiet 
subset (see § 3.2). Also, even the quasars in the high- 
redshift (z sw 4) faint sub-sample of Fan et al. (2001b) are 
quite massive (~ 10 9 M Q ). 

Kolmogorov-Smirnov tests confirm the apparent simi- 
larity seen in Figure 1 of the distributions in mass, lu- 
minosity, and Eddington ratio for all high-redshift sam- 
ples (§ 3.3), including the "color-selected sample" of Fan 




2 4 6 40 42 44 

Redshift log [L X (1350A or 5100A)/ergs s-'A" 1 ] 

Fig. 2. — Distributions of FWHM(H 0) (for BQS; open circles) and FWHM(Civ) (remaining quasars) with (a) redshift and (6) continuum 
luminosity, L A (5100A) (for BQS) and La(1350A) (for the z > 1 quasars). The horizontal, dashed line marks FWHM = 2000 kms" 1 , the 
commonly adopted cutoff for selecting quasars. Samples shown are: BQS (open circles); intermediate redshift radio-quiet (open triangles) 
and radio- loud (solid triangles) quasars; Constantin et al. sample (open squares); Anderson et al. sample (asterisks); Fan et al. sample and 
z « 6 quasars (solid points). 
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2 4 6 3.5 4 4.5 5 

Redshift Redshift 

Fig. 3. — Distribution of bolometric luminosities, Lboh as a function of redshift for (a) all the quasar samples considered here and (6) the 
z > 3.6 quasar samples only for clarity; error bars are omitted but the average error is shown in upper right corner. See Fig. 2 for symbols. 
The dashed curves represents the evolution of L? JLg) (lower) and L£ ol (10L^) (upper) with redshift. The evolution of L£ i at z > 2.5 is 
an extrapolation of the observed evolution at lower redshift (see text). The dotted curve segments denote limits imposed by sample selection 
and survey limits. The left panel shows limits for the z 2 quasars (V = 20 mag), the z 4 quasars (i* = 22.5 mag), and the z m 6 quasars 
(z* = 20.2 mag). The right panel shows the limit of the SDSS in general (i* = 22.5 mag) and, for reference, the limits of Fan et al. 's "bright" 
(i* = 20.2 mag) and "faint" (i* = 21 mag) samples. 



et al. (2001b). The low-redshift quasars of the BQS span 
a larger range in each of these parameters simply because 
less-luminous quasars are still bright and easily detectable 
at low redshift. 

Figure lc suggests that in every redshift range there is a 
significant number of quasars with Lboi/^Edd greater than 
unity. However, given the uncertainties in both Mbh and 
bolometric corrections, the reality of super-Eddington lu- 
minosities is unclear. Nevertheless, barring some large sys- 
tematic error, all the values of Lboi/^Edd are high, almost 
always greater than 0.1. This is not surprising given that 
we are sampling the brightest quasars in every redshift bin. 
Indeed, the lower end of the distribution in Lboi/^Edd is 
due to survey flux limits and the apparent upper envelope 
to the C IV line width for a given luminosity (Fig 2b). The 
upper end of the distribution is truncated by the lack of 
objects with very high luminosity and/or emission lines 
narrower than about 2000kms _1 . Since active nuclei are 
commonly identified by their broad lines, it is not clear 
whether the discovery surveys with their selection criteria 
(and crude spectral resolution in early surveys) would be 
able to identify quasars with such narrow lines; even the 
SDSS quasar identification is not well-defined at present 
(Richards et al. 2001). The upper cutoff in Lboi is however 
real (see e.g., § 4.2). 



4.2. Redshift Distributions 

Figures 3 and 4 show how Lboi and Mbh distribute with 
redshift. The curve segments are lower limits imposed by 
the survey limits or sample selection. In Figures 5 and 6 
the -Mbh, Lboh and ^boi/^Edd distributions with redshift 
are shown binned in different ways and with measurement 

6 The luminosity function determined by Pei (1995) is consistent with 



uncertainties omitted for clarity. 

The lower envelope of the observed Lboi values is 
easily explained by the survey or selection flux limits 
(Fig. 3). The bolometric luminosity corresponding to L* B 
and lOL^and their evolutionary tracks [i.e., L^ ol (L* B ,z) 
and L£ ol (10 L B , z)] are shown for reference, where L* B is 
the characteristic (break) luminosity of the quasar lumi- 
nosity function, $(Lb,z). For simplicity, the parameteri- 
zation of Boyle et al. (2000) is adopted here, namely: 

where Lb is the B broad-band luminosity, and a and (3 
are the slopes of the faint and bright end of the luminosity 
function, respectively. The L£ ol (z) values were computed 
as 10 16 x L B (z). The bolometric correction factor to Lb 
was determined from the updated quasar SED described in 
§ 2. The L B and its evolution 6 was determined for ^6000 
quasars by Boyle et al. (2000) between redshifts zero and 
2.3 from the LBQS and 2dF Quasar Redshift Survey data. 
Boyle et al. extrapolate L B (z) to z > 2.3 based on the fit 
to the data at lower redshifts. Figure 3 shows that the 
3.5 < z < 5 quasars studied here are very luminous: most 
are well above several L* B (see e.g., Fig. 3b), assuming 
L bol (L|j,z) can be trusted for short extrapolations (from 
z w 2.3). The extrapolation to z w 6 is much more un- 
certain, but, evidently, the most distant quasars detected 
are unusually bright. Also, the BQS quasars are clearly 
among the most luminous quasars at low redshift, as is 
well known, while the z ~ 2 quasars are, on average, L* B 
quasars or slightly more luminous. A preliminary analy- 
sis of the Large Bright Quasar Survey (LBQS) data from 
Forster et al. (2001) of quasars at 1.3 Z z S 3.0 suggests 
that the brightest of these quasars may exceed Lboi ~ 10 

that of Boyle et al. (2000) within the parameter uncertainties. 
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Redshift Redshift 

Fig. 4. — Distributions of Mbh with redshift for (a) all samples considered here, and (b) the 3.5 < z < 5.0 samples only for clarity; the 
average error is shown in upper right corner. See Fig. 3 for symbols. The dashed curves show the observed (2 < 2.5) and predicted (2 > 2.5) 
evolution of Afg H determined by assuming L^ ol = L^. dd (lower) and Lt .= 0.1 £ Edd (upper), respectively. The line segments denote limits 
in Mbh imposed by the survey flux limits (see Fig. 3) assuming FWHM(Civ) = 2000 km s -1 (1000 kms" 1 ) for z « 2 (z >4) quasars, and 
adopting the average (dashed segment) and minimum (solid segment) spectral slopes for the quasars at each redshift interval. 



ergss -1 and A/bh~ 10 10 M q . If indeed so, the relatively 
lower Lboi values of the intermediate-z sample (see Fig. 3) 
are consistent with the fact that these quasars were in- 
tentionally selected not to be the brightest at a given z 
(Vestergaard 2000). 

One possible explanation for the high luminosities is 
gravitational lensing by foreground galaxies too faint to 
be detected at present. The lensing fraction is strongly 
dependent on the slope (/?) of the bright end of the op- 
tical luminosity function at z > 3 (e.g., Wyithe & Loeb 
2002). Fan et al. (2001a) determined the luminosity func- 
tion for the well-defined "color-selected sample" of bright 
quasars at z w 4 and find a slope of (3 = 2.58, somewhat 
flatter than the slope, [3 — 3.43, of the bright end of the 
z < 2.3 luminosity function (Boyle et al. 2000). This shows 
that the z ss 4 "color-selected sample" contains a larger 
fraction of the most luminous quasars, in accordance with 
the findings here. More importantly, this bright end slope 
allows a handle on the lensing fraction at high-z, assum- 
ing the "color-selected sample" is representative of bright 
z > 4 quasars. Wyithe & Loeb (2002) find a relatively high 
probability (P) that the luminosity of z w 6 quasars is en- 
hanced by lensing: 7% < P < 30% for 2.58 < (3 < 3.43, 
respectively. The luminosities of the z m 4 quasars are 
not expected to be significantly lensed: 4% < P < 13%. 
These estimates show that very few of the high-z quasars 
studied here are likely to be gravitationally lensed. Even 
if as many as 13% of the most luminous z ss 4 quasars are 
lensed, the average Lboi of the remaining quasars is still 
well above the predicted L£ ol (z « 4): <L bo i>« 10 46 8 
ergss -1 = 9.4 x < L^ ol (z) >. Evidently, with the cur- 
rent survey flux limits, these high-z quasar samples only 



probe the brightest tail of the quasars possibly present at 
z w 4. Future surveys need to go much deeper 7 to prop- 
erly test whether or not the z > 3 luminosity function is 
significantly different than that at lower redshifts. 

Figure 4 shows that even the high-z black-hole masses 
are quite large 8 and, in particular, do not show a decrease 
relative to the masses of the intermediate-z quasar sam- 
ple. However, the masses of the z « 2 LBQS quasars may 
reach 10 10 M Q , as noted earlier. Thus, a slight mass de- 
crease may exist at z > 3.5; the apparent upper envelope 
in the quasar masses at 4 < z < 5 furthermore suggests 
this. Also shown in Figure 4 is the relative position of 
Afg H (z), determined from the evolution of L£ | assum- 
ing L£ ol = Lg dd oc Mg H . Since it is unclear whether all 
quasars radiate at the Eddington luminosity, the relative 
location of Mg H (z), computed assuming L£ ol = 0.1Lg dd , 
is also shown for reference. If Lboi = Lg dd ' s va hd, most 
of the quasars across the observed redshift range are very 
massive, reaching Mbh > 10 — 15 Mg H . This is true even 
at the lowest redshifts (z < 0.6). Regardless of the true 
value of Afg H (z w 4), the z > 3.6 quasars stand out as be- 
ing extremely massive. For most quasars, Lboi = 0.1 L Edd 
may be reasonable, given previous estimates at lower red- 
shifts (e.g., Wandel et al. 1999; Awaki et al. 2001; Czerny 
et al. 2001; but see also Collin et al. 2002). In this case, 
the BQS masses distribute almost evenly about the pre- 
dicted value of Mg H , with a range of a factor ~10 in each 
direction. Also, the values of Mbh f° r the z « 2 quasars 
are mostly in the range 0.1 — 1 x Mg H . Furthermore, 
if the A/g H (z) extrapolation to z w 4 is reasonable, then 
the z > 3.6 quasars have masses in the range ~ 0.1 — 10 
Mb H , similar to the BQS. The average observed Edding- 



7 An extrapolation of L* B shows that 1.0 (0.01) L* quasars are expected to have broad-band magnitudes: i* f«25 (30) mag at z ss 5 and z* ft: 
29 (31.5) mag at z ft 6. 

8 The mass estimates of the highest redshift (2=6. 41) quasar SDSS J114816.64+525150.3 based on both Mgn (Willott et al. 2003a) and Civ 
line widths (Barth et al. 2003) are in agreement (but note that individual mass estimates based on scaling relationships may be uncertain by 
a factor 10 (Paper I); the factor 3 uncertainty is valid for statistical samples only). 
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Fig. 5. — Distribution of Mbh with redshift shown binned in (a) Lbol an d (b) ibol/^Edd- Error bars are omitted for clarity. The labeled 
open star marks the likely position of an "average" broad-lined Lyman-break galaxy AGN (see § 5.2). 



ton ratio is, however, well above 0.1 (Fig. 1) — possibly 
closer to a value of 1 for the highest redshift quasars - 
with a typical Lboi/^Edd value of 0.4 — 0.5 at z > 1.5. 
Adopting this ratio would make Mg H fall about halfway 
between the two A/g H (z) curves shown in Figure 4, mak- 
ing the z » 2 quasars Mg H on average and the z > 3.6 
quasars as massive as 10 — 15 A/g H or even larger. Note 
that the absence of much less massive quasars, especially 
at high-z, is mainly due to the selection and survey lim- 
its (Figure 4). Also, the high Mbh values are not due 
to severe overestimates caused by the adopted method to 
estimate A/bh (§ 5.1). 

The very detection of the z > 3.6 quasars means they 
are very luminous, as already discussed. But, must they 
also necessarily be massive? Rejecting the possibility that 
they are grossly super-Eddington, the obvious possibili- 
ties for artificially enhancing their A/bh estimates are (a) 
gravitational lensing, and (b) relativistic beaming. Grav- 
itational lensing can account for some, but certainly not 
all quasars (30% at the very most), as argued above. The 
high-z quasars are also not likely to be strongly beamed, 
judging by the similar appearance of their UV spectra 
to those of low-z AGNs: beamed sources (BLLac and 
OVVs) have very weak or absent emission lines (e.g., Pe- 
terson 1997). The equivalent width (EW) of the UV 
lines in the high-z quasars are not much different from 
those of low-z AGNs. The average EW (Ly a +Nv) and 
EW(Ctv) in the Constantin et al. sample are 80A and 
~40A, respectively, while Fan et al. (2001b) find a mean 
£W(Lya+Nv) of 69 A for the "color-selected sample". 
These compare very well with the typical EWs in low- 
z AGNs: £W(Lya+Nv) = 75A and EW(Civ) = 35A 
(Peterson 1997). In fact, the EWs at low and high-z are 
remarkably similar considering more luminous quasars are 
expected to have lower line EW (the Baldwin effect; e.g., 
Baldwin 1977). So to conclude, since the z > 3.6 quasars 
are rather luminous, they must be massive as well. Ad- 
ditionally, the large emission-line widths (Figure 2) also 
indicate the black-hole masses are large. 



Figure 5 shows that even lower-luminosity (~ 
10 44 ergss _1 ) quasars can have relatively high masses, 
> 10 8 M Q (the BQS at z < 0.6). The tendency for 
-kboi/^Edd to increase with decreasing Mbh is particularly 
prominent for the z < 0.6 objects (Fig. 5b and 6b). No- 
tably, even the most massive z > 1.5 quasars also have 
high Eddington ratios (Fig. 6b). Also, a large fraction of 
the ATbh~ 10 8 — 10 9 M Q quasars appear to have super- 
Eddington luminosities (Figs. 5b and 6b). For the distant 
(z > 1.5) quasars, there is considerable overlap in luminos- 
ity for the two highest mass bins (Fig. 6a) and visa-versa 
(Fig. 5a). 

4.3. Overview of Results 

A number of conclusions can immediately be made 
based on the distributions described above. 

(1) At low redshift (z < 0.6), the brightest quasars, 
represented by the BQS quasars, exhibit a clear upper en- 
velope in Mbh (Fig. 4 and 5) and -L bo i (Fig- 3 and 6) 
which increases with redshift. Although such a sharp up- 
per cutoff is not seen for the more distant quasars there 
does appear to be a real physical ceiling of Mbh ~ 1O 1O M0 
and Lboi ~ 10 48 ergss -1 at all redshifts. The Lboi/^Edd 
ratio evenly covers the same approximate range (0.1 < 
-^boi/^Edd ~ 3) at all redshifts (Figs. 1 and 6). The lower 
envelopes of A/bh, ^boi, and iboi/^Edd for all the quasars 
are due to the original discovery surveys being flux limited. 

(2) The high-z quasars have relatively high values of 
A/bh and Lboi- The central masses of the highest redshift 
(z > 3.6) quasars do not drop significantly compared to 
the intermediate redshift (z w 2) objects. However, the 
z rs 2 quasars were not selected to be the brightest, and 
therefore the most massive, at this redshift (Figs. 1, 3, 
and 4). A preliminary analysis of the z > 1.3 LBQS sam- 
ple suggests that a mild decrease in, at least, Mbh may 
exist at z > 3.6, as the brightest z « 2 LBQS quasars may 
exceed Mbh~ 1O 1O M and L^oi ~ 10 47 ergss -1 . It is par- 
ticularly striking that the most distant quasars at z ~ 6 
have similarly high Mbh ~ 10 9 M Q as at lower redshifts. 
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Fig. 6. — Distributions of (a) Lbol an d (6) ^bol/^Bdd with redshift shown binned in Mbh- Error bars are omitted for clarity. The labeled 
open star marks the likely position of an "average" broad-lined Lyman-break galaxy AGN (see § 5.2). 



These quasars do have a high probability (P <30%) of be- 
ing gravitationally lensed (Wyithe & Loeb 2002). But if 
they are not lensed, the known z ~ 6 quasars must be very 
massive, ft has been argued that such massive quasars 
must necessarily reside in massive dark matter halos and 
be associated with very rare, high density peaks in the 
dark matter distribution at that epoch (e.g., Efstathiou & 
Rees 1988; Fan et al. 2001c; Volonteri et al. 2002). 

(3) Most of the objects have masses above Mg H (z), even 
at the lowest redshifts where L*(z) is reliable. Given the 
mean Lboi/i-Edd at z > 1.5, it follows that the actual A/g H 
in Fig. 4 probably lies between the two dashed curves. This 
makes the quasar masses of the intermediate- z sample on 
average Mg H and makes most of the z > 3.6 quasars 10 
— 15 Mg H , if the luminosity function can be reasonably 
extrapolated to z « 4. 

5. DISCUSSION 

It was argued in section 2 that black-hole masses of dis- 
tant AGNs and quasars can be predicted using their broad- 
line widths, their luminosities, and scaling relationships. 
For a statistical sample of objects this method appears ac- 
curate to within a factor of about 3. When this method 
is applied to the quasar samples described in § 3 it is con- 
cluded that quasars at 4 < z < 6 have black-hole masses of 
10 s - 10 10 M©, with an average of about 10 9 M Q (§ 4). In 
spite the fact that the currently known high-z quasars rep- 
resent only the upper end of the luminosity function (e.g., 
Fan et al. 2001a) (and presumably the upper end of the 
mass function), it is important to appreciate that black- 
hole masses as large as a few times 10 9 M Q are capable of 
forming by z = 6. Certainly, the high luminosities of these 
quasars suggest this must be the case (Fan et al. 2001c), 
but the present results provide evidence of high black-hole 
masses independently of arguments based only on the Ed- 
dington limit. 

9 Bechtold et al. (2003) overestimate the oqx values by adopting a UV slope of —0.3. Fan et al. (2001b) find that an average slope of —0.8 
is more appropriate for z ?s 4 quasars. For example, for two of the objects studied by Bechtold et al. the UV slopes are measured by Fan 
et al. and the resulting arjx are 1-47 and 1.48, much closer to the average of 1.43 measured at lower redshifts (Elvis et al. 2002; note that the 



It is currently unknown how these massive black 
holes form or on what time scale. As argued by Fan 
et al. (2001c), it takes 20(iboi/^Edd) _1 e-folding times to 
grow a 10M© black hole with radiative efficiency of 100% 
to a mass of ~ 5 x 10 9 M Q . As this time scale is within, 
although close to, the age of the Universe at z « 6 for a 
cosmology of Hq = 65 km s -1 Mpc _1 , Slo = 0.35, and A 
= 0.65, this growth scenario is not ruled out by its time 
scale. 

These conclusions depend on whether or not the scal- 
ing relations, defined for nearby AGNs and quasars, can 
be applied to luminous and very distant quasars. The 
main justifying reason is that the currently available data 
at radio, UV, and X-ray wavelengths show no clear evi- 
dence that the properties of z > 4 quasars are any different 
than those of other luminous quasars known at lower red- 
shift. First, quasar rest-frame UV spectra are very similar 
at all redshifts up to at least z w 4.5 (e.g., Constantin 
et al. 2002; Dietrich et al. 2002a; barring luminosity ef- 
fects such as the Baldwin effect). Also, as discussed in § 4, 
the values and distributions of A/bh, -^boii an d -^boi/^Edd 
are also similar, especially at z > 1.5. Second, at radio 
energies, the SDSS 3.6 < z < 5 quasars exhibit similar 
properties (~10% fractional occurrence and similar radio 
spectra) as known for z < 3 quasars (Carilli et al. 2001a; 
Stern et al. 2000; see Petric et al. 2003 for corroborating 
results for z > 5 quasars). Moreover, the X-ray properties 
of high- z quasars (e.g., the optical-UV to X-ray slope, aox, 
and X-ray slope, T x ) are probably not much different from 
those of their lower-z cousins. Brandt et al. (2002) and 
Mathur, Wilkes, & Ghosh (2002) independently concluded 
this for the aox values of three z ~ 6 quasars. Bechtold 
et al. (2003), on the contrary, argue that z > 4 quasars 
have relatively steeper aox values (> 1.5 with an average 
of 1.6±0.1) based on Chandra data of 17 optically selected, 
radio-quiet (3.7 < z < 6.3) quasars, but their aox values 
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are slightly overestimated 9 . Also, the aox slope extends to 
flatter values (aox ~ 1-2) when z « 4 quasars from other 
studies are included (Fig. 1 by Bechtold et al. ), similar to 
average aox values (1.55 < < aox > ~ 1-65) for z < 3 
radio-quiet quasars (Yuan et al. 1998). However, Vignali 
et al. (2003) also report steeper aox values for their z w 4 
quasars, but argue it is a luminosity effect. Nevertheless, 
Vignali et al. find the X-ray and broad-band spectra sim- 
ilar to those of low-z quasars. Bechtold et al. also deter- 
mine X-ray slopes of T x « 2.5 (1.5) for their 4.0 < z < 4.5 
(z > 5) quasars. These slopes are somewhat flatter than 
the steepest slopes of 3 to 4 observed for z < 2 ROSAT 
quasars. However, this apparent flattening of T x with red- 
shift is most likely a selection effect because different en- 
ergy ranges of the X-ray spectrum are observed at differ- 
ent redshift, as also suggested by Bechtold et al. (2003). 
In conclusion, the apparent similar properties across the 
SED of distant quasars (z > 1.5) indicate that the prop- 
erties of the central engine of the high-z luminous quasars 
are sufficiently similar that the scaling relationships can 
used to estimate black-hole masses of these objects. The 
efficacy of this method is discussed in more detail in § 5.1. 

The second key result presented in § 4, is the clear pres- 
ence of a ceiling to the Lboi and Mbh values at all redshifts. 
From redshift zero to 0.6, the maximum values of all three 
parameters increase to a certain value and then level off at 
higher redshifts. The upper envelopes in Mbh and Lboi are 
genuine, since none of the sample selections impose upper 
limits in the source brightness. Also, there are no apparent 
selection biases preventing a detection of very broad-lined, 
bright quasars, with the possible exception of the z « 2 
quasars (§ 4). The apparent lack of very luminous and very 
massive quasars may be due to either or both of the follow- 
ing explanations. First, very luminous (massive) quasars 
are expected to be very rare, given the steep decline of 
the upper end of the luminosity (mass) function. Second, 
the above-mentioned ceilings may be due to some phys- 
ical limit to black hole growth. Possible causes for such 
a maximum sustainable mass and luminosity include: (a) 
limits in the dark matter gravitational potential, (b) limits 
in the fuel supply rate, (c) transition of accretion mode or 
timescale, and (d) changing accretion disk (i.e., emission) 
properties with black-hole growth (e.g., Caditz, Petrosian, 
& Wandel 1991; Yi 1996; Kauffmann & Haehnelt 2000). 

Two recent studies (Bechtold et al. 2003; Netzer 2003) 
report quasar masses at high-z (z « 4 and z i 3, respec- 
tively) reaching or exceeding 10 10 M Q in apparent conflict 
with the above conclusion. However, these Mbh values 
are not inconsistent with the current analysis once differ- 
ences in cosmology and statistical uncertainties are taken 
into account, as explained in the following. Firstly, Net- 
zer (2003) finds a maximum mass of 10.4 dex assuming an 
R-L slope of 0.7 and H = 70 km s^Mpc" 1 , fl m = 0.3, 
and f^A = 0.7. For this cosmology, his Mbh values are 
about a factor 2 (or < 0.35 dex) higher than for the cos- 
mology adopted here (Ho = 75 km s _1 Mpc _1 , q = 0.5, 
and A = 0) . Combining this with the probable uncertainty 
in his mass values of —0.6 dex, given a measurement error 



of < 0.25 dex and a statistical uncertainty in the method 
of —0.5 dex (Paper I), the highest mass values of Netzer's 
sample and those analyzed here are statistically consistent. 
An important point to keep in mind is that individual Mbh 
values can be in error by as much as a factor 10. The 
fraction (< 9%) of quasars in Netzer's sample with Mbh 
> 10 10 M Q is also within the expected fraction (—10%) 
of objects for which the Mbh values may be off by 1 dex 
or more (Paper I). Secondly, the maximum Mbh value 
(-10.7 dex) reported by Bechtold et al. (2003) for H = 
50 km s _1 Mpc _1 , q = 0.5, and A = 0, corresponds to 
— 10.45 dex in the current cosmology. Although, Bechtold 
et al. do not quote their errors, the uncertainties in their 
Mbh values are presumably as large as or even larger 10 
than those of Netzer (2003), and, hence, their maximum 
Mbh value is similarly consistent with the upper boundary 
< 10 10 M Q , seen in the current study, to within the errors. 

Netzer (2003) also claims that M BH £ 1O 1O M is incon- 
sistent with the locally defined Mbh —<t relationship (a 
is the stellar velocity dispersion in the host galaxy bulge) , 
since the implied a — 700km s" 1 is not observed locally. 
As Netzer notes, either (1) the Mbh values are overesti- 
mated, because, for example, the R — L relationship does 
not extrapolate to high-z quasars, or (2) the Mbh — o re- 
lationship changes with redshift. It is argued in § 5.1 that 
the R — L relationship is not likely to be misleading and 
that the scaling relationships used here yield reasonable 
Mbh estimates to within a factor <4. Furthermore, the 
lack of local black-holes with masses ~ 10 10 M Q does not 
mean that the large Mbh values inferred from the scal- 
ing relationships are strongly overestimated. The reason 
is that in addition to a possible evolution of the Mbh — c 
relationship and the statistical uncertainty in the Mbh es- 
timate, the lack of local black-holes with masses ~ 1O 1O M 
can possibly be explained by the relatively small volume 
probed in the local Universe: for the comoving space den- 
sity of luminous, massive (Z 5 x 1O 9 M0) quasars at z w 4 
of < 12.5 Gpc~ 3 (Vestergaard & Osmer, in preparation), a 
volume approximately 20 times larger than that currently 
probed locally (— 100 Mpc) is needed to detect one such 
massive quiescent black-hole. 

5.1. Reliability of the Mass Estimates for High Redshift 

Quasars 

It is a serious concern whether the scaling relationship, 
defined for a representative sample of nearby AGNs and 
quasars, can rightfully be applied to luminous, distant 
quasars. Specifically, can extension of this relationship 
cause the masses of the distant quasars (~ 10 9 M Q ) to be 
severely overestimated? Netzer (2003) also points out that 
the high masses he finds at z < 3 appear inconsistent with 
the local Mbh — & relationship. In this section, the efficacy 
of the scaling relationships for luminous, distant quasars 
is discussed. It is argued that application of the scaling 
relationships to such quasars is valid, owing mainly to the 
similarity of quasar spectra across the known redshift and 
luminosity ranges. In particular, the scaling relation is 

radio-quiet quasars observed with ROSAT have < qqx >= 1.55; Yuan et al. 1998). It may thus not be the norm that high-rcdshift quasars 
have very steep optical-X-ray slopes; but see Vignali et al. 2003). 

10 This is owing to their crude correction of their FWHM(Clv) values to FWHM(H/3), so to apply the optical scaling relationship of Kaspi 
et al. (2000) to obtain central masses, and their adoption of rather flat UV slopes (—0.3) relative to the steep slopes (—0.8) measured by Fan 
et al. (2001b) for their representative z ss 4 quasars (see footnote above). 



Massive Black Holes at High Redshift. 



11 



unlikely to be inaccurate by a factor >4 in the mass. 

The mass estimates will obviously be misleading if the 
assumptions made when extending the scaling relation- 
ships beyond the nearby AGNs break down. The two as- 
sumptions for M BH oc FWHM(Civ) 2 AL°- 7 (1350A) are: 
(1) that the size, i?BLR, of the (Civ) emitting region in 
the broad line region - i.e., its distance from the contin- 
uum source - can be estimated reliably from the contin- 
uum luminosity, L\(1350A), and (2) that the line width 
is a reasonable proxy for the gravity-dominated velocity 
dispersion of the line-emitting gas. The critical issues are 
therefore: (a) whether or not the R — L relationship is 
valid at higher redshift and for more luminous quasars 
than the sample for which it was defined, and (b) whether 
or not the C iv emission gas is virializcd or, for example, 
the Civ profile has strong blue asymmetries which may 
be a sign of BLR outflows as suggested for the narrow-line 
Seyfert 1 (NLS1) galaxies (e.g., Leighly 2000). Each issue 
is addressed in turn below. 

First, the R — L relationship is not likely to be much 
different for the more luminous and more distant AGNs. 
The reasons are twofold: (1) the R — L relationship is 
not extrapolated very far in luminosity: for the highest 
ALa(1350A) values of 47.4 dex, the R — L relationship is 
extrapolated by less than 1.5 dex, a small fraction of the 
luminosity range of >4dex, over which the R — L rela- 
tionship is defined (10 42 < Ai A (5100A)/ergss" 1 < 10 46 ), 
and (2) AGN and quasar spectra look very similar at all 
rcdshifts and luminosities considered here, i.e., line flux ra- 
tios, line equivalent widths, and SEDs are all similar (e.g., 
Elvis et al. 1994; Kuhn et al. 2001; Constantin ct al. 2002; 
Dietrich et al. 2002a, 2002b; see also § 2). Specifically, if 
the R — L relationship does not apply to luminous, distant 
quasars, the spectra of these objects would look somewhat 
different than they do. This can be quantified by photoion- 
ization modeling. In essence, if the scaling relationship 
(eq. 4) results in a mass overestimate of, for example, a 
factor 10 for a quasar with Lboi ~ 10 47 ergss _1 , the virial 
theorem shows that the Ly a and C iv emission lines are 
emitted so close (<33 light-days) to the ionizing source that 
these lines arc produced very inefficiently, and, in particu- 
lar, C in] emission cannot be generated. But, the presence 
of C in] in all spectra of quasars (at least up to z « 4.3; 
Dietrich et al. 2002a) confirms the existence of low-density 
BLR gas subject to a lower ionizing flux from which Lya 
and Civ are also much more efficiently emitted 11 : this 
gas will necessarily be located at a greater distance from 
the ionizing source, close to the distance inferred from the 
luminosity. 

Second, there is no evidence that the kinematics of the 
BLR in high-z AGNs are not dominated by gravity. This 
conclusion is based on the following 3 arguments. 

1. Most, if not all, of the BLR gas 12 is in virial mo- 
tion around the central source, including the high- 
ionization UV lines: all the AGNs with suitable 
multiple emission line measurements demonstrate a 
robust virial relationship between their line widths 

11 These conclusions are based on the line emission resulting from gas with a range of column densities subject to the ionizing flux from a 
source with Lbol ~ 10 47 crgss _1 as quantified by the photoionization grids of Korista ct al. (1997). 

12 This has so far been established for the following emission lines: SilvA1400, ClvA1549, HenA1640, Cm] A1908, H/3A4861, HellA4686 (see 
e.g., Peterson & Wandel 1999). 

13 Some spectra are noisy and the reality of such a prominent triangular line shape is difficult to discern. 



and variability time lags (Peterson & Wandel 1999, 
2000; Onken & Peterson 2002); while only four 
AGNs are testable so far, the general similarity of 
AGN spectra (i.e., line shape, line widths, and vari- 
ability properties) argues that this relationship ex- 
tends to most AGNs. 

2. The distribution of Civ line widths in spectra of 
distant, luminous quasars is similar to that for 
nearby AGNs (typical values ~4500kms _1 with a 
range from 2000 - 10000km s" 1 in FWHM), includ- 
ing those for which the scaling relations are cal- 
ibrated. Also, the Civ widths of the nearby BQS 
quasars are consistent with their H (3 line widths, as- 
suming virial BLR kinematics and that H /3 is emit- 
ted twice as far from the central source as C iv, as 
established for three of the four AGNs discussed 
in item 1 (see also Korista et al. 1995; see com- 
ment below on the one exception). Furthermore, 
the similar FWHM(Civ) distributions, where most 
are narrower than ^8000kms _1 (e.g., Figure 2), 
argue that the unknown, but possible, contribu- 
tion to single-epoch line widths from non-variable 
optically thin emission (Shields, Ferland, & Peter- 
son 1995; Peterson 1997) is highly unlikely to make 
FWHM(Civ) in error by a factor >2; such an er- 
ror translates to a mass uncertainty of a factor <4, 
comparable to the statistical uncertainty (§ 2). 

3. Outflow of the BLR line-emitting gas is not im- 
portant, at least for the quasars considered here. 
This is supported by two facts: (i) none of the 
Civ profiles of the z > 1.5 quasars (§ 3.2, § 3.3), 
and most quasars in general, clearly 13 resemble 
the broad, triangular, blue asymmetric C iv pro- 
files of the high luminosity NLSls; those Civ pro- 
files have been argued to signify strongly outflow- 
ing high- ionization line-emitting gas (e.g., Leighly 
2000), as NLSls are very luminous for their Mbh, 
and (ii) while a small fraction (<20%) of the 
Civ profiles in the high-z sample (§ 3.3) show 
weak, blue-ward asymmetries (mostly suspected 
due to broad, blended Niv]A1486), those quasars 
do not have much larger masses and luminosities 
than the remaining quasars: a conservative es- 
timate is an 0.2 dex increase on average, domi- 
nated by the lower-quality spectra of the Ander- 
son et al. sample; the 'asymmetric subset' of the 
Constantin et al. and Fan et al. quasars exhibit av- 
erage deviations of ± 0.1 dex. Also, the C iv profiles 
of the z w 2 quasars are not typically asymmetric, 
and any asymmetries are restricted to the very base 
of the profile and do not affect FWHM. Moreover, 
the strongest asymmetries do not extend blueward 
(Vestergaard 2000; Vestergaard, Wilkes, & Barthel 
2000). 

In conclusion, blue asymmetries and UV outflows are gen- 
erally not a concern for the Mbh estimates of quasars. 
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Nonetheless, extreme care should be exercised when us- 
ing the C iv line width to estimate the mass when signs 
of significant outflows are present. For this reason, mass 
estimates based on Civ are not well suited for the most 
luminous NLSls. 

It is important to emphasize that although all AGNs 
so far tested show the BLR is virialized (i.e., the virial 
product Rv 2 is constant for all emission lines measured), 
for one object (3C 390.3) Civ is narrower than H/3 and 
has a larger lag in contrast to the other three objects. 
This behavior is not understood, but 3C 390.3 may have 
an unusual BLR structure as it is the only broad-line radio- 
galaxy with monitoring data and it, moreover, has double- 
peaked Balmer lines (Peterson & Wandel 2000). Neverthe- 
less, it is very important for the general applicability of the 
scaling relations that future studies address this issue. 

5.2. Do Black Holes "Mature" Before Their Host 
Galaxies Do? 

Given that black holes with large masses appear as early 
as z « 6 and given the intimate connection between the 
formation of galaxy bulges and their central, supermas- 
sive black holes inferred from the Mbh —o~ relationship, 
the question naturally arises whether or not this close con- 
nection also exists at high rcdshift. Specifically, do these 
"mature" black holes (i.e., those that have reached, say, 
~ 10 8 — 10 9 M Q ) also reside in "mature" host galaxies (i.e., 
virialized systems with the majority of their stars already 
formed) with aged (>lGyr) stellar populations? Among 
the currently available data on high- z AGN host galaxies, 
there is considerable circumstantial evidence that massive 
active black holes may, in fact, reside in galaxies that are 
not fully relaxed, virialized systems and where many of the 
stars are yet to be formed. This is based on (a) the in- 
ferred high star-formation rates (SFRs) from the existing 
large masses of cool dust, (b) the presence of significant 
amounts of molecular gas in some sources, (c) the presence 
of large-scale young starbursts in at least one well-studied 
z w 4 radio-galaxy with a large dust mass, (d) the strik- 
ingly different morphology of host galaxies of high- z radio- 
quiet quasars and radio galaxies compared to their lower- z 
cousins, and (e) the presence of massive, active black-holes 
in small, young galaxies at z > 3. Each issue is discussed 
in turn below. 

5.2.1. The Evidence 

Massive Star Formation in High-z Quasars High- 
z AGN host galaxies are apparently experiencing massive, 
large-scale star formation: the far-IR/sub-mm SEDs of 
z Ri 4 quasars (e.g., Benford et al. 1999; Carilli, Menten, 
& Yun 1999; Carilli et al. 2001a) are similar to those of 
ultra-luminous IR galaxies (ULIRGs) and starburst galax- 
ies (e.g., Downes & Solomon 1998; Genzel et al. 1998; 
Rowan-Robinson 2000; Klaas et al. 2001). ULIRGs are 
dominated by extreme rates of star formation (>1OOOM0 
yr ) that yield unique far-IR/sub-mm SEDs, which are 
independent of whether or not a central AGN is present 
(e.g., Downes & Solomon 1998; Klaas et al. 2001). The 
rather cool (T<100K; typically - 40 - 50 K) starburst- 
heated dust dominates the far-IR emission (e.g., Rowan- 

14 The inferred dust masses are, in fact, lower limits as cold dust (T 
there is probably a factor ~3 more gas (e.g., Dunne et al. 2000). 



Robinson 2000; Dunne et al. 2000; Klaas et al. 2001), but 
AGN-heated dust (T >100K) emits in the near-IR, short- 
ward of — 60^m (e.g., Wilkes et al. 1999b; van Bemmel & 
Dullemond 2003; Figure 8a by Farrah et al. 2003). Large 
amounts (> 10 8 M©) of cool (<~ 40 — 50 K) dust are com- 
monly inferred 14 from the rcstframc far-IR/sub-mm emis- 
sion from z w 4 quasars (e.g., Benford et al. 1999; Car- 
illi et al. 2001a; Omont et al. 2001; Priddey & McMahon 
2001; Isaak et al. 2002). Large dust masses are often taken 
to imply high SFRs, -1000 - 2000 M Q yr" 1 (Benford 
et al. 1999; Carilli et al. 1999, 2000, 2001b, 2003); such 
high SFRs are even inferred for the undetected sources 
(Isaak et al. 2002). Also, non-thermal radio emission con- 
tributes little to the shape of the cm to sub-mm SED, 
confirming its origin in massive star formation (Carilli 
et al. 2000, 2001a, 2001b; Yun et al. 2000). 

The far-IR properties of the luminous z « 4 SDSS 
quasars studied here appear very similar to those of other 
z w 4 quasars with far-IR data and to far-IR selected 
AGNs. Most of them have radio to far-IR SEDs consis- 
tent with starburst heating of large dust masses (Carilli 
et al. 2001a), and the inferred dust masses (Md us t) and 
SFRs are similar to other z « 4 quasars: Figure 7 shows 
that for the subset of 30 quasars with Mbh and L\, Q \ es- 
timates, the 13 quasars detected at 1.2 mm all have Mbh 
> 10 9 Af Q . Even for the undetected quasars, the upper 
limits cannot rule out the existence of SFRs of order 300 
- 700 M Q yr" 1 . 

This line of argument is not without caveats. For 
example, it can be argued (e.g., Sanders et al. 1989; 
Kuraszkiewicz et al. 2003) that the IR SEDs of nearby 
quasars can be explained purely by AGN heating, although 
this explanation is not universally accepted (e.g., Farrah 
et al. 2003; Klaas et al. 2001). Nevertheless, no correla- 
tion is observed between quasar optical-UV luminosity 
and far-IR luminosity, as otherwise expected if AGN heat- 
ing dominates the cooler dust (e.g., McMahon et al. 1999; 
Omont et al. 2001; Isaak et al. 2002; Priddey et al. 2003). 
Also, Archibald et al. (2002) argue that a high dust mass 
may indicate a galaxy is in a late starforming phase (>75% 
complete), as discussed below. 

Molecular Gas. Huge amounts of cold molecular gas 
(> 10 11 M Q ) are present in z > 3 far-IR luminous sources. 
This is (a) determined from CO emission line measure- 
ments (Ohta et al. 1996; Omont et al. 1996; Guilloteau 
et al. 1999; Carilli et al. 1999, 2002a, 2002b; Papadopou- 
los et al. 2000), and (b) inferred from the dust masses and 
typical dust-to-gas ratios obtained from low-z starbursts 
(e.g., Sanders, Scoville, & Soifer 1991; Dunne et al. 2000). 
Such vast gas reservoirs show that large amounts of stars 
are yet to be formed (e.g., Solomon et al. 1997; Downes & 
Solomon 1998), and are evidence for ongoing massive star 
formation, especially when combined with large masses of 
cool dust (e.g., Ohta et al. 1996; Ivison et al. 1998; Frayer 
et al. 1999; Papadopoulos et al. 2000; Carilli et al. 2002a, 
2002c). In some cases, the cold gas clearly constitutes 
a significant fraction of the dynamical mass in the sys- 
tem, emphasizing the rather young nature of the galaxy 
(e.g., Ohta et al. 1996; Guilloteau et al. 1999; Solomon 

K) is a poor emitter (e.g., Omont et al. 2001; Chapman et al. 2003); 
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et al. 1997; Downes & Solomon 1998; Tacconi et al. 1999). 
For some high-z sources, the sub-mm emission is (a) clearly 
separated (tens of kpc) from the AGN nucleus, and/or 
is (b) distinctly extended (e.g., Papadopoulos et al. 2000; 
Omont et al. 2001; Isaak et al. 2002; Carilli et al. 2001a, 
2002b, 2003). This is strongly suggestive of galactic- 
scale starforming regions and renders an AGN origin less 
likely (Carilli et al. 2003). Furthermore, when CO obser- 
vations of high-z sources are conclusive, the high masses of 
cold gas inferred from dust mass estimates are confirmed 
(Hughes et al. 1993; Ohta et al. 1996; Omont et al. 1996; 
Papadopoulos et al. 2000). Non-detections of CO line 
emission, however, do not necessarily rule out a young 
starburst interpretation because the observations are still 
not very sensitive 15 . 



Young Stellar Populations. For one far-IR/sub-mm 
luminous high-z AGN, the z w 3.8 radio galaxy 4C 41.17, 
stellar populations synthesis modeling has established the 
presence of large-scale young starbursts with ages in the 
range 0.01 Gyr to 0.4 Gyr with a most probable age of 
0.07 Gyr (Dey et al. 1997; Dey 1999). This is also impor- 
tant because it confirms the existence of starburst activity, 
expected based on the strong far-IR emission alone. Al- 
though it is difficult to generalize based on one object, it 
is pertinent to emphasize that such analysis is very diffi- 
cult to perform for most high-z AGNs, given the strong 
nuclear glare in quasars and the faintness of galaxies at 
high-z. Nevertheless, it is important for future studies to 
establish whether or not high-z far-IR luminous AGNs in 
general have young starburst activity. 

Archibald et al. (2002) infer starburst ages of 0.5 — 
0.8 Gyr for 4C 41.17 from a monochromatic luminosity ra- 
tio and their models of dust production and depletion with 
time in a young galaxy, and conclude the galaxy is in a late 
starforming phase. These ages are clearly at odds with 
the results of the more sophisticated age determination 
discussed above. Further developments in this area seem 
necessary to resolve these apparent inconsistencies. More 
importantly, the measured dust mass of 4C 41.17 is also 
consistent with a starburst age in the range 0.075 Gyr to 
0.095 Gyr (Fig. 1 by Archibald et al. 2002), in good agree- 
ment with the results of Dey (1999). In conclusion, while 
dust contents alone is an inaccurate measure of starburst 
ages, large dust masses are not uniquely tied to >lGyr 
old galaxies but also appear in galaxies with starbursts as 
young as 0.04 Gyr (for M dust > 10 8 M Q ; Fig. 1 by Archibald 
et al. 2002). 

It is therefore an intriguing possibility that some high-z 
AGN host galaxies may have rather young stellar pop- 
ulations. However, a comparison with the restframe IR 
emission is necessary to determine the extent to which 
these young stars dominate the total stellar mass in place 
at that epoch. The observed apparent extension of the 
K — z relation, defined for z < 1.6 radio galaxies (Lily & 
Longair 1984), to z > 3 radio galaxies and quasars is some- 
times used to argue for somewhat old stellar populations 
in their host galaxies. However, it is important to appreci- 
ate that the increasing scatter in the observed relationship 

15 The success-rate of CO detection is currently relatively low due to 
band widths compared to the sometimes broad line widths (Ohta ct al 
gas is too cold, the CO line transitions are very weak. 



at z > 2 is also consistent with galaxies comprising mainly 
young stellar populations forming at 4 < z < 5 (De Breuck 
et al. 2002): the stellar mass in the host galaxies (quanti- 
fied by the X-band magnitude) is too low in many cases to 
place all high-z radio-galaxies on the K — z relation with 
insignificant scatter. Also, it has been argued that the na- 
ture of proto-galaxies or young, less massive galaxies may 
also conspire to fall along the K — z relation at z > 3 (Eales 
& Rawlings 1993; Willott et al. 2003b). Hence, it seems at 
least plausible that a good fraction of high-z host galaxies 
may be relatively young, even if this interpretation is not 
universally accepted. 

Host Galaxy Morphology While quasar host galaxies 
appear fully assembled by z « 1, AGN host galaxy proper- 
ties become increasingly more difficult to quantify at z > 1 
(Kukula et al. 2001). Practically, radio galaxies are the 
only AGNs we can currently study beyond z « 3, owing to 
the lack of nuclear glare that affects optical observations of 
quasar hosts. Radio-galaxy hosts do not appear fully as- 
sembled and dynamically relaxed. Specifically, in contrast 
to their lower-z cousins with mostly single, smooth com- 
ponents, z > 3 radio galaxies display very clumpy struc- 
ture: multiple, small star-forming components are embed- 
ded in large-scale (^50 kpc) faint emission, indicative of 
galactic systems in early stages of assembly (e.g., Graham 
et al. 1994; van Breugel et al. 1998; Pentericci et al. 1999, 
2001; Papadopoulos et al. 2000). Optically obscured z > 2 
radio galaxies also appear to be in an early stage of for- 
mation (Reuland et al. 2003). 

For high-z radio-quiet quasar hosts, the best-quality 
imaging data available are likely those from NICMOS 
HST studies of 2 < z < 3 sources (Ridgway et al. 2001; 
Kukula ct al. 2001). In contrast to radio-galaxy hosts, 
these are relatively structureless, single components. Also, 
at 2 < z < 3, radio-quiet hosts are often more compact 
(scale-lengths <4kpc and ^2.3 kpc on average; Ridgway 
et al. 2001) than (a) radio-loud hosts at that epoch (scale- 
lengths ~ 11 — 18 kpc; e.g., Best et al. 1998; Kukula ct 
al. 2001; Ridgway et al. 2001), and than (b) quasar hosts 
at z w 0.2 (mean scale-lengths of ^8.2 kpc; Bahcall et al. 
1997; McLure et al. 1999). Moreover, the hosts of z « 2 
radio-quiet quasars are apparently rather young: Ridgway 
et al. find the 2 < z < 3 (Mb ~ —24 mag) radio-quiet 
quasar hosts contain only ~10% — 20% of their final stel- 
lar mass. A caveat is that the Ridgway et al. results may 
not reflect the properties of the 2 — 3 magnitude more lu- 
minous radio-quiet quasars observed at z « 4, but deep 
imaging of z > 3 radio-quiet quasar hosts is currently un- 
available. 

Can AGNs Occur in Young Galaxies? Probably the 
most convincing evidence that relatively young galaxies at 
z > 3 are capable of hosting massive, active black holes 
comes, somewhat surprisingly, from the distant, small- 
scaled, young, star-forming Lyman-break galaxies (LBGs; 
e.g., Steidel et al. 1996, 1999). The presence of an ac- 
tive black hole can be inferred for a tiny fraction (~1%) 
of the ^1000 LBGs detected in a survey at z « 3 (Stei- 

sensitivity limits of the detectors and the relatively narrow accessible 
. 1996; Papadopoulos et al. 2000; Omont 2003). Also, if the molecular 
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Fig. 7. — Distribution of inferred star formation rate (SFR) 
and mass of dust with estimated black-hole mass for a sub- 
set of the 2 « 4 SDSS quasars analyzed here. The far-IR 
emission emitted by the cool dust was observed by Carilli 
ct al. (2001a). Detections at A b s = 1.2 mm are shown as 
solid circles, while the open triangles denote 3 a upper lim- 
its. The dust masses and SFRs are here computed using 
the prescription of Omont et al. (2001). Specifically, SFR 
~ 360 <5mf 5gp (5250 GHz /mJu) Mq yr" 1 , where 5gp is the 
ratio between the starburst's total luminosity and the far- 
IR luminosity, and <5mf is "a function of the present mass 
composition of the stellar population" (see Omont ct al. for 
explanation) and has values in the range 0.8 < 5mf ~ 2. 
Since starburst heating dominates the far-IR, <5sf > 70% is 
assumed, and thus 0.6 S <5mf<5sf ~ 1-4; a fair estimate is 
^MF^SF ~ 1 yielding a lower limit on the SFR. 



del et al. 2002) by the presence of broad emission lines in 
their rest-UV spectra. Specifically, the average spectrum 
of these broad-lined LBGs (Steidel et al. 2002) suggests 
that the mass of the black hole is typically rather high 
(> 10 8 M Q ): Table 1 lists estimates of the average values 
of Mbh, iboi, and -Lboi/^Edd for different cosmologies and 
assumed dust extinction factors. Based on the uncorrected 
observed fluxes, the broad-lined LBGs have masses Mbh 
« 1O 8 M , lower than the higher luminosity quasars at any 
rcdshift beyond 1.5 (Figure 5), and they are about an or- 
der of magnitude less luminous (Figure 6). In spite of the 
large uncertainties, the inferred average Mbh value is in- 
triguingly large considering (a) the relatively faint nature 
of both the nuclear and host galaxy emission, (b) the star- 
formation activity in the host galaxy (Steidel et al. 1996), 
and (c) the low stellar mass in the galaxies at z w 3: they 
typically only contain ^10% — 30% of their expected final 
stellar mass (Ridgway et al. 2001; Papovich et al. 2001). 
This suggests that the massive central black hole in the 
LBGs formed early and fast compared to the underlying 
stellar component. 

A possible caveat is the uncertainty as to whether or 
not these faint LBG-type AGNs have similar SEDs to 
the bright quasars, on which both the Mbh and, espe- 
cially, Lboi estimates are based. But among the lower- 
luminosity AGNs in the nearby Universe, only the faint, 
very low-luminosity LINERs display dramatically differ- 
ent SEDs (Ho 1999), mainly missing the "big blue bump" 
emission from the accretion disk. Scyfert galaxies, with 
which the LBGs with broad high-ionization lines may be 
most similar, display only minor SED differences in the 
optical-UV region compared to the more luminous quasars 
(e.g., Sanders et al. 1989; Elvis et al. 1994). Moreover, the 
relatively large Civ line widths (4670±500kms _1 ; very 
similar to the mean FWHM(C iv) of the luminous z « 4 
quasars) and line equivalent widths (see Figure 1 by Stei- 
del et al. ) support a high central mass (as argued in § 4.2), 
in spite the faint nature of the galaxies. Therefore, the de- 
rived mass and luminosity values are crude, but probably 



not too unrealistic. 

5.2.2. Discussion and Summary 

Although largely circumstantial, the evidence presented 
above suggests that the z > 3 massive, active black holes 
typically reside in relatively young (possibly dynamically 
unrelaxed) galaxies in their major star-forming phase. The 
implication is thus that the black hole either formed long 
before the stars in the galaxy or sometime during the 
first star-formation phase and thereafter grew so fast rel- 
ative to the stellar mass that the nucleus shined long be- 
fore the galaxy was fully formed. Rix et al. (2001) and 
Omont et al. (2001) independently reached similar conclu- 
sions based on different arguments. Rix et al. find z ~ 2 
radio-quiet quasars to have either higher nuclear luminosi- 
ties or fainter host galaxies than lower- z AGNs. They ar- 
gue that the most plausible explanation is that the ratio 
of Mbh to the stellar mass increases with redshift. Omont 
et al. (2001) use millimeter data on z « 4 quasars to ar- 
gue for a higher black-hole growth rate relative to the 
surrounding bulge stellar component in the host galaxy 
at early epochs, consistent with the current study and 
that of Rix ct al. Such an evolutionary scenario is consis- 
tent with some hierarchical models (e.g., Haehnelt & Rees 
1993). Also, Silk & Rees (1998) suggest that early black- 
hole growth might in fact help regulate star formation in 
the host galaxy and its evolution. The overall conclusion 
reached here is thus that black holes likely reach maturity 
before the stellar populations in their host galaxies. 

Archibald et al. (2002) claim that quasar activity is de- 
layed 0.5 — 0.7 Gyr from the onset of star formation, 
during which >75% of the stars are formed, in apparent 
contradiction to the results presented above. However, 
Archibald et al. assume that black holes grow by accretion 
with a radiative efficiency (n) of 10%, while the value of n 
and its constancy with time (or redshift) are highly uncer- 
tain. A lower radiative efficiency allows faster black-hole 
growth, and indeed, some recent estimates imply r\ « 0.03 
(Schirbcr & Bullock 2003). Furthermore, one can equally 
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Table 1 

Estimated mass, luminosity, and Eddington ratio of the average 
broad-lined agn among lyman-break galaxies. 



Extinction 
Factor 
Applied 
(1) 


Ho 
(kms x ) 
Mpc -1 ) 

(2) 


logL„(1350A) 
(erg s _1 Hz -1 ) 

(3) 


log M B h 

(M ) 

(4) 


log £boi 
(erg s- 1 ) 

(5) 


-^bol/iEdd 

(6) 


None 


50 


29.54±0.02 


8.2±0.1 


45.56±0.2 


0.20±0.10 


None 


75 


29.19±0.02 


7.9±0.1 


45.20±0.2 


0.15±0.08 


50 


50 


31.24±0.02 


9.4±0.1 


47.26±0.2 


0.64±0.32 


50 


75 


30.89±0.02 


9.1±0.1 


46.90±0.2 


0.50±0.25 



a These estimates and measurement uncertainties are based on the composite spectrum of Steidel 
et al. (2002). While it is based on non-spectrophotometric data, it should after all be representative 
of the mean brightness level of these AGNs (C. Steidel, 2002, private communication), sufficient for the 
discussion here. Specifically, the typical C IV line width is measured (from the published plot and similarly 
to the other data analyzed here; § 3) to be 4670±500 kms -1 , and the 1350A flux is (2. Of ±0.09) ^Jy. Since 
the composite spectrum shows evidence of some dust reddening, the Mbh, £boi, and Lboi/iEdd values were 
also computed after making a correction for dust extinction. Steidel et al. (2002) estimate an extinction 
factor of ~ 50 for their narrow-lined AGNs, but do not quote such a factor for the broad-lined AGNs 
among their LBGs. Therefore, the estimate for the narrow-lined AGNs were adopted here in attempts to 
bracket reality somewhat. In addition to the listed values of Ho, qo = 0.5 is assumed. 



well expect black-hole formation and early mass growth to 
be very fast, as expected for the first generation of mas- 
sive stars (e.g., Abel, Bryan, Norman 2002). Silk & Rccs 
(1998) find it unlikely that the first black holes formed 
with M B h< 10 6 M Q ; in that case, it would take <0.2Gyr 
to generate a 10 8 M Q black hole from a 10 6 M@ seed ac- 
creting with r\ <10%, consistent with quasar activity ex- 
isting in very young galaxies, as it takes 0.3 Gyr to form 
50% of the stars in the galaxy (Archibald et al. 2002). 

It is pertinent to re-emphasize that the presence of 
young starbursts in high-z quasars cannot be inferred from 
their high dust content alone: a young evolving galaxy 
maintains a large dust mass (> 10 8 M Q ) from 0.04 Gyr 
to 1.5 Gyr after the initial starburst (e.g., Archibald 
et al. 2002). But the combination of high masses of cool 
dust and cold (molecular) gas is much more convincing. 
Therefore, it is important for future investigations to study 
the cold gas contents in high-z AGN host galaxies to test 
whether or not they are in the process of building up the 
majority of their stars. Furthermore, better coverage of 
the radio through near-IR SEDs of high-z quasars is re- 
quired to constrain the relative contributions of AGN and 
starburst heated dust for more accurate estimates of SFRs. 

6. CONCLUSIONS 

The main results of this work can be summarized as 
follows: 

(1) It is confirmed that the high-z quasars have cen- 
tral masses of at least a few times 10 9 M Q as previously 
estimated from their luminosities (Fan et al. 2001c). The 
masses, bolomctric luminosities, and Eddington ratios of 
the z > 3.6 quasars are very similar to those of their lower- 
redshift cousins. Relative to the evolved Lboi and Mbh 
the z w 4 quasars are much more luminous and more mas- 
sive than the typical quasars at z < 3.5. This assumes 



that the luminosity function, based on the Large Bright 
Quasar Survey and the 2-Dcgree Field Quasar Redshift 
Survey, and its extension to z « 4 is reasonably reliable 
(§ 4.2). 

(2) The scaling laws from which the black-hole mass is 
estimated are fully applicable to distant luminous quasars. 
This is argued based on the similarity of quasar and AGN 
spectral properties (line widths, line equivalent widths, 
line ratios, SEDs) at low and high redshift. In par- 
ticular, photoionization models confirm that if the ra- 
dius — luminosity relationship does not extend to high-z 
quasars, their spectra would look different than they do. 
Mass estimates, based on scaling relationships, are not 
likely inaccurate by a factor greater than ^3 to ^4 (§ 5.1) 
when applied to statistically significant samples. 

(3) There is a real ceiling to the Mbh and _L b oi values 
at all redshifts: M BH < 10 10 M Q and L bol < 10 48 ergss _1 . 
These upper limits are consistent with the recent stud- 
ies of Netzcr (2003) and Bechtold et al. (2003) to within 
the uncertainties and differences in adopted cosmologies. 
Moreover, the large M B h values at high-z are not inconsis- 
tent with local M B h values given the much smaller volume 
probed locally (§ 5). 

(4) Black holes reach maturity (i.e., reach masses Mbh 
« 10 8 - 1O 9 M ) before their host galaxies: AGN host 
galaxies at z > 3 do not typically appear fully assembled, 
dynamically relaxed, and/or the majority of their stars 
are not formed or are somewhat evolved (>1 Gyr). Since 
the evidence for relatively young quasar-host galaxies at 
z > 3 is mostly indirect, this interpretation needs to be 
tested further. For example, complete coverage of the ra- 
dio to IR SEDs of large representative samples of z > 3 
quasars will help provide improved constraints on the rel- 
ative starburst and AGN contributions to the far-IR emis- 
sion. Detailed studies confirming the correlation between 
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large dust masses, large reservoirs of cold molecular gas, 
and high star-formation rates in larger quasar samples are 
also required. No doubt, the combination of SIRTF and 
the planned FIRST/Herschel and ALMA telescopes will 
play an important part in shedding light on these issues. 
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